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We try to single out some qualitative new effects of the 
coupling to the (5-isovector-scalar meson introduced in a mini¬ 
mal way in a phenomenological hadronic field theory. Results 
for the equation of state {EOS) and the phase diagram of 
asymmetric nuclear matter {ANM) are discussed. We stress 
the consistency of the 5-coupling introduction in a relativistic 
approach. 

New contributions to the slope and curvature of the sym¬ 
metry energy and the neutron-proton effective mass splitting 
appear particularly interesting. 

A more repulsive EOS for neutron matter at high baryon 
densities is expected. 

Effects on new critical properties of warm ANM, mixing 
of mechanical and chemical instabilities and isospin distilla¬ 
tion, are also presented. The 5 influence is mostly on the 
isovectorlike collective response. 

The results are largely analytical and this makes the physi¬ 
cal meaning quite transparent. Implications for nuclear struc¬ 
ture properties of drip-line nuclei and for reaction dynamics 
with Radioactive Beams are finally pointed out. 


BAGS numbers: 21.65.-bf, 21.30.Fe, 25.70.Pq, 24.10.Jv 

I. INTRODUCTION 

Hadronic effective field theories [quantum hadrody- 
namics, QHD\ represent a significant improvement in 
the understanding of static and dynamical properties of 
nuclear matter and finite nuclei, described as strongly in¬ 
teracting systems of baryons and mesons. The main ap¬ 
plication of this approach has been the Relativistic Mean 
Field theory (RMF) p|J^, extremely successful in nuclear 
structure studies |^|-^ 

In recent years, due to the new possibilities opened by 
the Radioactive Beam facilities, the interest has moved 
towards a microscopic description of asymmetric nuclear 
systems, from the stability of drip-line nuclei to the phase 
diagram of asymmetric nuclear matter. The connection 
to astrophysics problems, supernovae explosions and neu¬ 
tron stars, is also quite evident il- The isovector 
channel has been introduced through a coupling to the 
charged vector p-meson. In the Hartree approximation 
of the RMF approach it leads to a simple linear increase 
of the symmetry term with the baryon density, without 
neutron/proton effective mass splitting |]l^. A complete 


study of the corresponding phase diagram for asymmet¬ 
ric nuclear matter has been performed in ref. 0 , with 
a thorough discussion of the new qualitative features of 
the liquid-gas phase transition. 

In this work we focus our attention on the introduction 
of the coupling to the isovector scalar channel, through 
the exchange of a virtual charged (5[ao(980)] meson. In 
this way, as well know in the isoscalar channel, we re¬ 
cover the genuine structure of the relativistic interac¬ 
tions, where one has the balance between scalar (attrac¬ 
tive) and vector (repulsive) ’’potentials”. Actually the 
(5-meson exchange is an essential ingredient of all NN 
realistic potentials and its inclusion in the QHD scheme 
has been already suggested p^, also on the basis of a 
relativistic Brueckner theory . Some first results 

of structure calculations for exotic nuclei are showing the 
importance of the (5-dynamics for the stability conditions 
of drip-line nuclei |T5|Jl^ ]. 

Aim of this paper is to present, within a RMF ap¬ 
proach, the main expected effects of the 5-field on sym¬ 
metry properties of the nuclear system, from EOS to 
n,p-mass splitting, and in particular on the nuclear re¬ 
sponse in unstable regions. We remind that the study 
of symmetric and asymmetric nuclear matter properties 
under extreme conditions is of great relevance for the 
understanding of the nuclear interaction in the medium. 
In particular a liquid-gas phase transition may occur in 
warm and dilute matter produced in heavy-ion collisions 
(multifragmentation events). This is one of the interest¬ 
ing open problems in theoretical and experimental nu¬ 
clear physics. Isospin effects on heavy ion collisions (frag¬ 
mentation and collective flows) have been largely dis¬ 
cussed in recent years, see the reviews in [^,Q. In this 
work we will try to offer a unified picture of the 5-meson 
exchange influence on nuclear structure and dynamics. 

In order to simplify the analysis and to pin down the 
most direct effects of the 5 contributions we will follow a 
quite reduced version of the RMF approach, yet includ¬ 
ing all the isoscalar (ct, w) and isovector (5, p) fields. Non¬ 
linear self-interaction terms are introduced only in the 
isoscalar scalar cr-channel, essential in order to get right 
incompressibility parameter at normal density po [0 , ^ ■ 
Such approximation scheme of minimal self-interacting 
terms is actually physically well justified in the baryon 
density range of interest here, roughly up to 2.5po [ 01 ■ 
Moreover it is consistent with the RMF assumptions of 
neglecting retardation and finite range effects in the field 
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dynamics. 

We like to stress the presence of several points of inter¬ 
est on the isovector channel results in this baryon density 
region [p^ : 

• The slope of the symmetry term, the Symmetry 
Pressure, just below pQ is of relevant importance for 
nuclear structure, being directly linked to the thick¬ 
ness of the neutron skin in n-rich (stable and/or un¬ 
stable) nuclei, see the recent discussions in [p5|-p7[ . 
Jointly to the n — p effective masses this informa¬ 
tion appears essential for the stability of drip-line 
nuclei 

Dissipative reaction mechanisms with asymmetric 
ions seem to be also quite sensitive to the same 
quantity p8|j2|]. 

• The nature of chemical instabilities and in general 
critical properties of warm and dilute asymmetric 
nuclear matter |]ll| , ^o| -p^ . 

• Isospin effects on the dynamics of heavy ion colli¬ 
sions at intermediate energies, see [Q. Reaction 
studies represent a very sensitive tool to test trans¬ 
port properties of symmetry effects, n,p chemical 
potentials and effective masses ||3^-|^. 

In this paper, we first determine the model parameters 
by fitting the properties of the symmetric and asymmet¬ 
ric nuclear matter at T = 0. We then extend the inves¬ 
tigation to finite temperature. In particular we derive 
the new boundary regions for mechanical and chemical 
instabilities and their dependence on the various mesons 
entering in the theory. 

This paper is organized as follows. In Sec. 2, the 
equation of state for nuclear matter at finite temperature 
is derived. Symmetry energies and effective masses are 
discussed in Sec.3, where the neutron matter EOS is 
also evaluated. The mechanical and chemical instabilities 
are studied in Sec. 4. General comments as well as a 
summary of the main conclusions are presented in the 
last Section. 


II. EQUATION OF STATE FOR NUCLEAR 
MATTER AT FINITE TEMPERATURE 

The starting point is the relativistic Lagrangian den¬ 
sity of an interacting many-particle system consisting 
of nucleons, isoscalar (scalar a, vector cu) and isovector 
(scalar S, vector p) mesons: 


£ = - (Mat - ga(j) - gsf ■ 5) 

- gpi^T 


+ 1^% ■ ■ d>^6- 

( 1 ) 

Minimal self-interacting terms, as discussed in the in¬ 
troduction, are included only in the cr-channel. 

Here (j) is the cr-meson field, the w-meson field, 
the charged p meson field, 6 the isovector scalar field of 
the J-meson. We define = 

dpbu — dvbp. U{(j)) is the nonlinear potential of the a 
meson : U{(j)) = \a(p + \bp. 

The field equations in mean field approximation 
(RMF) are: 

- {Mn - ga4 - gsTsSs) - gu>pi^a 
-gpPnbo)^ = 0 
ml(t) + acp + bcj)^ = 44 = gaPs 
= gu>4l°4 = gu:PB 

"m-lbo = gp4pT34 = gpPBs. 

mp3 = gs4T34 = gsps3, (2) 


where pB3 = PBp - PBn and ps3 = PSp - psn, PB and ps 
being the baryon and the scalar densities, respectively. 

Neglecting the derivatives of mesons fields, in the 
mean field approximation the energy-momentum tensor 
is given by 

= i4lp.9u4 + \=ml(lP- + U{4>) + 

- ]^mlbxb>']gp,,. (3) 

The properties of nuclear matter at finite temperature 
are described by the thermodynamic potential Q. ^From 
statistical mechanics for a system in a volume V we define 
Q = —pV = —^InZ H], where /3 the inverse of tempera¬ 
ture, /3=1/T, and Z is the grand partition function given 
by Z = Here H is the hamilto- 

nian operator, Bi and pi are respectively nucleon num¬ 
ber operators and thermodynamic chemical potentials, 
{i = p,n). The equation of state for nuclear matter at 
finite temperature can be obtained from the thermody¬ 
namic potential Q. 

The energy density has the form 

/ p\3p 1 

—^E*{k){n4k) + n^k)) + 

L—tL^p '' 

+U{4)+ ^mlujl + ^mlbl + ^mpl, ( 4 ) 

and the pressure is 


P = 


i—n^p 


d^k 


{2n4Et{k) 

1 


(ni(fc) -I- ni(fc)) 


-Ui4) 


2 2,^ 2 I 2 

+ -m„ba 




^ 2 i2 

ImjSl (5) 
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where Ei* = \J. The nucleon effective masses 
are defined as 


M* = Mjy — g„(j) =F (— proton, + neutron). ( 6 ) 


The ni{k) and ni(k) in Eqs. (4)-(5) are the fermion and 
antifermion distribution functions for protons (i = p) and 
neutrons {i = n): 

l + exp{(£;,*(fc)-/r*)/T} ’ 

and 

l + ex.p{{E*{k) + p*)/T}' 

The effective chemical potentials p* are given in terms 
of the vector meson mean fields 

Mi = T gpbo (- proton, + neutron), (9) 


where pi are the thermodynamical chemical potentials 
Pi = dejdpi. At zero temperature they reduce to the 

Fermi energies Epi = kp- + M*"^. 

The same results can be directly obtained from the ex¬ 
pectation value of the energy-momentum tensor, show¬ 
ing the thermodynamical consistency of the mean field 
approximation |^. 

By using the field equations for mesons, the equations 
of state for thermal matter can be rewritten as 




i—n,p 


A'ik ^ 

^Et{k){nm+n,{k)) + -mlcf^ 




9Z 


2nil 


Pb 


9p 2 
PB3 


2ml 


9l 2 

2MS3j 


2m 


( 10 ) 


and 


P = 


i—n,p 


d^k 


(27r)3 E;*(/c) 


(ni(fc) +ni{k)) - 


-U{cj>) 


2 m 2 


Pb + 


2 m^ 


Pb3 


^.2 

2PS3- 


2m 


( 11 ) 


We remind that the baryon densities ps are given by 
(7 is the spin/isospin multiplicity) 

/ USjL 

^(n(fc)-h(fc)), ( 12 ) 

while the scalar densities ps are 

/ c\^k M* 

^ — {n{k)+nik)). (13) 

At the temperatures of interest here the antybaryon 
contributions are actually negligible. 


In order to study the asymmetric nuclear matter, we 
introduce an asymmetry parameter a defined as a = 
(PBn - PBp)/PB = 

The energy density and pressure for symmetric and 

asymmetric nuclear matter at finite temperature and the 

n,p effective masses can be calculated self-consistently 

from Eqs. ( 6 ) to (13) just in terms of the four boson 
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coupling constants, fi = (^), i = o',ui,p,6 and the two 
parameters of the a self-interacting terms, A = -^ and 



The isoscalar meson parameters are fixed from sym¬ 
metric nuclear matter properties at T = 0: satura¬ 
tion density po = 0.16fm~^, binding energy E/A = 
— 16MeE, nucleon effective mass M* = 0.75Mn (Mjv = 
939MeV) and incompressibility Ky = 240MeE at po- 
The fitted f^, fui, A, B parameters are reported in Table 
I and have quite standard values for these minimal non¬ 
linear RME models. 


Table I. Parameter sets. 


parameter 

Set I 

Set II 

NL3 

fc ifm^) 

10.33 

same 

15.73 

U ifm^) 

5.42 

same 

10.53 

fp ifm'^) 

0.95 

3.15 

1.34 

fs ifm^) 

0.00 

2.50 

0.00 

A ifm-^) 

0.033 

same 

- 0.01 

B 

-0.0048 

same 

-0.003 


In the table we report also the NL3 parametrization, 
widely used in nuclear structure calculations Q. We re¬ 
mind that the A^L3-saturation properties for symmetric 
matter are chosen as po = 0.148 /to“3^ = Q.QMn, 

Ky = 271.8Mey and the symmetry parameter (see next 
section) is 04 = 37.AMeV. 

The corresponding phase diagram (pressure isotherms) 
is shown in Fig.l. The local minimum disappears at 
the critical temperature Tc, determined by ^|tc = 

= 0 - The obtained critical temperature of sym¬ 
metric matter is Tc = 15.86 MeV. 
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FIG. 1. Pressure as a function of the baryon density for 
symmetric nuclear matter (a=0) at different temperatures. 


III. THE SYMMETRY ENERGY 


T = 0. Then we discuss first the case at zero temper¬ 
ature, in order to fix the isovector coupling constants. 
^From Eqs.(10,15) we can easily get an explicit expres¬ 
sion for the symmetry energy 0 


E, 


sym 


I \ — ^ I 1 r _ 1 f _ 

~ 2^^ El\\ 


M* 


‘Pb 


■ fsA^kp, M*)] 
(16) 


where kp is the nucleon Fermi momentum correspond¬ 
ing to Pb, Ep = -I- M*^) and M* is the effective 

nucleon mass in symmetric NM^ M* = Mjy — g^ip. 

The integral 


A{kp,M*) = 


{2ny 


(Ek 


(fc2 +M*2)3/2 


(17) 


has a simple analytical structure which makes quite 
transparent the effect of the (5-meson on the symmetry 
energy. With some algebra we can get 


Symmetry properties of asymmetric nuclear matter 
are univocally fixed from the coupling constants of the 
isovector channels Experimentally we have just one 
relatively well known quantity, the bulk symmetry en¬ 
ergy 04 of the Weiszaecker mass formula, in the range 
30 — 35 MeV. In all the RMF models with only the 
isovector p meson this assigns the fp constant, that gives 
also the ’’slope” of the symmetry term (the Symmetry 
Pressure) [|I^, apart small non-linear contributions [^ . 

When we further include the (5-meson, from the 04 - 
value we will fix only a combination of the two coupling 
constants, while the S itself will imply interesting new 
contributions to the slope and curvature of the symme¬ 
try energy and a neutron-proton effective mass splitting. 
These will be the main points discussed in the following. 
We note that most results are analytical and this will 
largely improve the physical understanding of the effects 
due to the coupling to the S isovector scalar channel. 

The symmetry energy in asymmetric NM is defined 
from the expansion of the energy per nucleon E{pB,Oi) 
in terms of the asymmetry parameter 

E{pB,a) = ^ = E{pb) + Esym{pB)a^ 

Pb 

+0(a4) + ... (14) 


and so in general 

_ ld‘^E{pB,a) 


E, 


da^ 


_ 1 I 

|a —0 2^^ d ^ lpS 3—0 

' B3 


The Bulk Symmetry Parameter 

The bulk symmetry parameter is just the value of 
Esym corresponding to the normal density pB = Po, at 


A{kp,M*) = 3{^-^) (18) 

making use of 

PsiT = 0) = ^[kpEp - ln(^^^^)] 

Pb(T = 0) = ^4 (19) 

Expanding the scalar density in terms of (kp/M*)'^ [|j, 
the quantity Eq.(18) can be written as 


A{kp,M*) 


3pB A , kp 2 _3 fcf 4 _5 fcf g 


that can be used to derive a similar expansion for the 
symmetry energy eq.(16). 

We see that A{kp,M*) is certainly very small at low 
densities, below po- It can be still neglected up to a 
baryon density pB — 3/9o, where it reaches the value 
0.045/m“^ (with our symmetric NM parameters), i.e. 
a correction of about 10% in the denominator of Eq.(16). 

Then in the density range of interest here we can use, 
at the leading order, the much simpler form of the sym¬ 
metry energy, extremely nice in order to get the main 
qualitative new features due to the 5-meson coupling: 


JJ, I X 1 kp I !(/• r \2l 

EsymiPB) 2 ~ /< 5 ( ) \PB 


( 20 ) 


We see that, when the S is included, the observed 04 value 
actually assigns the combination [fp — fs{^)^] of the 
{p, 6) coupling constants. If ^ 0 we have to increase 
the p-coupling (see Fig.l of [|2|). In our calculations 
we use the value 04 = 30.5MeV. In Table I the Set I 
corresponds to fs = 0. In the Set II fs is chosen as 
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2.5/m^. Although this value is relatively well justified, 
we stress that aim of this work is just to show the 
main qualitative new effects of the i5-coupling. 

In order to have the same 04 we must increase the p- 
coupling constant of a factor three, up to fp = 3.15/m^. 
From the structure of Eqs.(16, 20) it is clear that there 
is a connection between the scalar field (isoscalar, cr, and 
isovector, 5) contributions, since both are acting on ef¬ 
fective masses. A strong cr-coupling (smaller M*, e.g. 
see NL3 parametrization) can compensate a strong 5- 
coupling. In our evaluation we keep fixed the cr param¬ 
eters {fa,A,B) leading to M* = 0.75Mn. In any case 
a larger value of Egym at high baryon densities, due to 
the relativistic mechanism discussed in the following, see 
Fig.2, will be always present when the <5 is included. 

At subnuclear densities, pB < Po, in both cases, (p) 
and (p -I- ^), from Eq.(20) we have an almost linear de¬ 
pendence of Egym on the baryon density, since M* Ep 
as a good approximation . Around and above po we ex¬ 
pect a steeper increase in the (p -I- <5) case since M*/Ep 
is decreasing, see Fig.2. 



Pb/ Po 

FIG. 2. Total {kinetic + potential) symmetry energy as 
a function of the baryon density. Dashed lines: (p). Solid 
lines: {p + 5). Upper curves: zero temperature. Lower curves: 
T = 8MeV. 

This is an interesting aspect to look at more in detail 
since it actually represents a general relativistic effect 
due to the coupling to scalar mesons. In some sense it is 
the equivalent, in the isovector channel (5 vs. p), i.e. for 
the symmetry energy, of the saturation mechanism we 
have in the isoscalar channel (tr vs. uj) for the symmetric 
matter. 

The scalar charged 5-meson, like the neutral tr, is act¬ 
ing on the nucleon effective masses, introducing a n,p 
splitting, see Eq.( 6 ). This causes a negative contribution 
to Egym, Eqs.(16),(21), since it reduces the gap between 
n, p Fermi energies due to the different Fermi momenta in 
asymmetric NM. In fact in a n-rich matter the neutron 


Fermi momentum increases while the neutron effective 
mass decreases, see Eq.( 6 ) and Fig .6 (the opposite for 
protons). 

Such negative 5 contribution is reduced at high densi¬ 
ties due to the ’’Lorentz contraction” factor {M*/Ep)'^, 
Eqs.(16) and (20), that in general gives the attenuation 
of the scalar interactions with increasing baryon density. 

We stress the consistent picture of a symmetry energy 
built from the balance of scalar (attractive) and vector 
(repulsive) contributions, with the scalar channel becom¬ 
ing weaker with increasing baryon density. This is in¬ 
deed the isovector counterpart of the saturation mecha¬ 
nism occurring in the isoscalar channel for the symmet¬ 
ric nuclear matter. From such scheme we get a further 
strong fundamental support for the introduction of the 
5-coupling in the symmetry energy evaluation. 



Pb/Po 

FIG. 3. p- (open circles) and 5- (crosses) contributions 
to the symmetry energy, second and third terms of Eq.(16). 
The dashed line is the approximate 5-contribution of Eq.(20), 
see text. The solid line is the linear extrapolation of the low 
density behaviour, plotted to guide the eye. 

In Fig.3, using our parametrizations, we quantitatively 
show the interplay of the two contributions (p and 5) to 
Egym, second and third terms of Eqs.(16),(20), as a func¬ 
tion of the baryon density. The crosses (negative 5 con¬ 
tribution) follow a linear behaviour up to roughly po and 
then they tend to saturate due to the Lorentz contrac¬ 
tion factor. In correspondence we will see the increase of 
the total Egym shown in Fig.2. We note the accuracy of 
the approximate 5 contribution given by Eq.(20), dashed 
line in Fig.3. 

In Fig.4 we show the Equation of State (energy per 
nucleon) for pure neutron matter {a = 1 ) obtained with 
the two parameter Sets, I (p) and II (p -I- 5). 

The values are in good agreement with recent non- 
relativistic Quantum-Monte-Carlo variational calcula¬ 
tions with realistic 2— and 3—body forces 0 - The in¬ 
clusion of the 5-coupling leads to a larger repulsion at 
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baryon densities roughly above l.Spo- This could be of 
interest for the structure of neutron stars and the possi¬ 
bility of a transition to new forms of nuclear matter. 



FIG. 4. EOS for pure neutron matter. Dashed line: (p). 
Solid line: {p + 6). 


is actually fixed by the 04 value) and a genuine S ex¬ 
tra contribution of about 20MeV from the second term. 
When we include also the kinetic part (from the first term 
of Eq.(20)) we have a total slope parameter going from 
L(p) = +84:MeV to L(p -h (5) = +108MeV. 

We note again that the slope parameter, or equiva¬ 
lently the Symmetry Pressure Psym = PoL/i, is of great 
importance for structure properties, being linked to the 
thickness of the neutron skin in n-rich (stable and/or un¬ 
stable) nuclei , and to the assessment of the drip¬ 

line I p^ . Moreover the same parameter gives an estimate 
of the shift of the saturation density with asymmetry (at 
the lowest order in a^) 


,^From the previous analysis we expect to see also in¬ 
teresting i5-effects on the slope (symmetry pressure) and 
curvature (symmetry incompressibility) of the symmetry 
energy around poj of relevant physical meaning. This will 
be the main subject of the following discussion. 


Symmetry Pressure and Symmetry 
Incompressibility 


In order to have a quantitative evaluation of the effect 
on the slope and curvature, and also to compare with 
other estimates, we use the expansion of the symmetry 
energy around po 

r, / N L pg—pQ KsymfPB~P0\2 , 

Esym(PB) =a4 + -( - -H —r|-(- + ■■■ 

3 Po 18 Po 

( 21 ) 


with L and Ksym respectively related to slope and cur¬ 
vature of the symmetry energy at po 


l = I„. 


K, 


dpB 

-^svm 


and 


= 9p^( 


dpi 


( 22 ) 


that can be easily obtained from a linear expansion 
around the symmetric value po(q; = 0 ) j^ ]. 

Eq.(24) has a simple physical meaning: in order to 
compensate the symmetry pressure in asymmetric mat¬ 
ter we have to move the zero of the P(pb) curve to lower 
densities. The amount of the shift will be inversely pro¬ 
portional to the slope, given by the incompressibility. 



Pb/ Po 

FIG. 5. Pressure as a function of the baryon 

density at zero temperature for various asymmetries 
a = 0 (bottom),0.5 (middle),0.8 (top). Dashed lines: (p). 
Solid lines: (p + 5). 


,^From Eq.(20) we get a potential contribution to the den¬ 
sity variation of Egym given by (after some algebra) 


dEgym I _ Iff f P] \ 
dpB 


Ep 


1 PB dM* 
3 ~ Jd* dpB 


(23) 


Around normal density po the first term is fixed by the 
04 value and so we get a net increase of the slope, due 
to the S, given by the second term (always positive since 
(dM*/dpB) < 0 ). 

,^From the first term of Eq.(23), with our parametriza- 
tion, we get a potential contribution to L of 45MeE (this 


All that can be seen in Fig.5 where we report the total 
pressure P(pb) for various asymmetries, at zero temper¬ 
ature, in the two cases, p-meson only (dashed lines) and 
(p -h 5) (solid lines). 

It is instructive to perform a similar discussion for the 
curvature parameter Kgym, Eq.(21). Now the potential 
contribution is exclusively given by the 5-meson, with a 
definite positive sign, as we can see from the previous 
discussion. This will be a large effect on the total since 
the kinetic part is quite small ||4^ . If we compare with 
non-relativistic effective parametrizations | p4|j4^ ], when 
we add the 5-meson we move from a linear to a roughly 
parabolic ps-dependence of the symmetry energy. 
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With our parameters we pass from a Ksym{p) = 
+7MeV (only kinetic) to a Ksym{p + (5) = +12QMeV. 
So this quantity appears extremely interesting to look at 
experimentally, as recently suggested from reaction mea¬ 
surements 1^^. The problem is that the effect on the 
total incompressibility of asymmetric matter, that likely 
could be easier to measure, is not trivial. We can eval¬ 
uate a shift of the incompressibility with asymmetry, at 
the same o? order as in Eq.(24), as 

AKv{a) = {Ksym-6L)a^ (25) 

So what really matters for the total incompressibility is 
the combination {Ksym — 6T), with the possibility of a 
compensation between the two terms. Just by chance this 
is actually what is happening in our calculations since for 
the above combination we get —A97MeV in the case of 
only p coupling, and —SOdMeE when we add also the 5. 
Indeed we see from Fig.5 that the slopes of the pressure 
around equilibrium density (P(po) = 0) are very close in 
the two cases, (p) and {p + 5). 


Effective Mass Splitting 

Another qualitative new result of the J-meson coupling 
is the n, p-effective mass splitting in asymmetric matter 
Ip^ , see Eq.(6). In Fig.6 we report the baryon density 
dependence of the n,p effective masses for a = 0.5, {N = 
iZ) asymmetry, calculated with our Set II parameters, 
compared to the symmetric case 



Pb/ Po 

FIG. 6. Neutron, proton effective masses vs. the baryon 
density for a = 0.5 {N = 2>Z). The dashed line corresponds 
to symmetric nuclear matter. 

We get a splitting of the order of 15% at normal density 
POj and increasing with baryon density. Unfortunately 
from the present nuclear data we have a very small knowl¬ 
edge of this effect, due to the low asymmetries available. 


This issue will be quite relevant in the study of drip¬ 
line nuclei. Moreover we can expect important effects 
on transport properties ( fast particle emission, collec¬ 
tive flows) of the dense and asymmetric NM that will be 
reached in Radioactive Beam collisions at intermediate 
energies. 

The sign itself of the splitting would be very instruc¬ 
tive. As we can see from Eq.(6) in n-rich systems we 
expect a neutron effective mass always smaller than the 
proton one. The same is predicted from more micro¬ 
scopic relativistic Dirac-Brueckner calculations |^. At 
variance, non-relativistic Brueckner-Hatree-Fock calcula¬ 
tions are leading to opposite conclusions § 13 . Still in 
the non-relativistic picture quite contradictory results are 
obtained with the Skyrme effective forces. The most re¬ 
cent parametrizations, SLy-type of Skyrme forces 
give the proton effective mass above the neutron one, in 
agreement with our calculations. Previous Skyrme-like 
forces, instead, yield a splitting in the opposite direction, 
but also show unpleasant behaviours in the spin channel 
(collapse of polarized neutron matter, see discussion in 
@)- 

We like to note again that, as shown before, a decreas¬ 
ing neutron effective mass in n-rich matter is behind the 
relativistic mechanism for the symmetry energy, balance 
of scalar (attractive) and vector (repulsive) contributions 
in the isovector channel. 


Finite Temperature Effects 

In the temperature range of interest in this paper, be¬ 
low the critical temperature of the liquid-gas phase 
transition, of the order of 15 — 16 MeV (see Fig.l), 
temperature effects on the symmetry properties are not 
expected to be large. Indeed the contributions of an¬ 
tifermions, that could modify all the terms with the 
scalar densities, see Eq.(13), are still very reduced. In 
Fig.2 we present also the symmetry energies calculated 
at r = 8MeV, lower solid and dashed curves. For both 
cases, p and {p + 6) the variation is quite small. We have 
a reduction mainly coming from the kinetic contribution 
due to the smoothing of the n,p Fermi distributions. We 
note that this result is in full agreement with relativistic 
Brueckner-Hartree-Fock calculations . 

The effect is even smaller for the mass splitting, given 
by a difference of scalar densities. In the next section we 
will study in detail the phase diagram of heated asym¬ 
metric nuclear matter, focussing in particular on the in¬ 
stability regions. 
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IV. MECHANICAL AND CHEMICAL 
INSTABILITIES 


dfip 

dy 


< 0 , 


(28) 


Heavy-ion collisions can provide the possibility of 
studying equilibrated nuclear matter far away from nor¬ 
mal conditions, i.e. to sample new regions of the NM 
phase diagram. In particular the process of multifrag¬ 
mentation allows to probe dilute nuclear matter at finite 
temperatures, where in the symmetric case, see Fig.l, we 
expect to see a phase transition of first order of liquid-gas 
type, as suggested from the very first equations of state 
built with effective interactions 

The multifragmentation phenomenon may be inter¬ 
preted as a signal of such a phase transition in a fi¬ 
nite system, when the nuclear matter in the expansion 
phase enters the region of spinodal (mechanical) instabil¬ 
ity. Nuclear matter is however a two-component system 
consisting of neutrons and protons: binary systems have 
more complicated phase diagrams due to the new concen¬ 
tration degree of freedom. In particular for asymmetric 
nuclear matter a qualitative new feature in the liquid-gas 
phase transition is expected, the onset of a coupling to 
chemical instabilities (component separation) that will 
show up in a novel nature of the unstable modes, the 
mixture of density and charge fluctuations leading to an 
Isospin Distillation effect |pd] , 011 ■ Indeed equilibrium 
thermodynamics, as well as non-equilibrium kinetic con¬ 
siderations predict that an asymmetric system will sepa¬ 
rate into more symmetric larger fragments (” liquid”) and 
into neutron-rich light fragments (”gas”). So the chem¬ 
ical instability can be investigated experimentally just 
by measuring the N/Z ratio, or isospin content, of the 
fragments. 

Since the effect is driven by the isospin dependent part 
of the nuclear equation of state, here we will look at the 
influence of the S coupling on this new liquid-gas phase 
transition. 

It is known that the stability condition of a two- 
component, n — p, thermodynamical system is given by 


dP 

dp 


T,v 


f dp.p \ 

\ 9y J rp p 


> 0 , 


(26) 


where P is the pressure, Pp is the proton chemical poten¬ 
tial and y the proton fraction Zj related to the asym¬ 
metry parameter a = 1 — 2y. 

Eq.(26) is equivalent to set the free energy to be a 
convex function in the space of the n,p density oscil¬ 
lations, SpnjSpp. In charge symmetric matter isoscalar 
(total density) Spn + Spp and isovector (concentration) 
6pn — Spp oscillations are not coupled and we have two 
separate conditions for instability 


mechanical^ i.e. 
region), and 



(27) 


vs. density oscillations (the spinodal 


T,P 


chemical, i.e. vs. concentration oscillations. 

In asymmetric matter the isoscalar and isovector 
modes are coupled and the two separate inequalities do 
not maintain anymore a physical meaning, in the sense 
that they not select the nature of the instability. Inside 
the general condition Eq.(26) the corresponding unstable 
modes are a mixing of density and concentration oscilla¬ 
tions, very sensitive to the charge dependent part of the 
nuclear interaction in the various instability regions . 

In dilute asymmetric NM (n-rich) the normal unsta¬ 
ble modes for all realistic effective interactions are still 
isoscalarlike, i.e. in phase n — p oscillations but with 
a larger proton component This leads to a more 

symmetric high density (liquid) phase everywhere under 
the instability line defined by Eq.(26) and consequentely 
to a more neutron-rich gas {Isospin Distillation). Such 
’’chemical effect” is driven by the increasing symmetry 
repulsion going from low density to roughly the satura¬ 
tion value and so it appears quite sensitive to the symme¬ 
try energy of the used effective interaction at subnuclear 
densities. It could provide a good opportunity to differ¬ 
entiate the various EOS isospin dependences |0| . 

In this section we study the effect of the 5-coupling on 
the instability region given by Eq.(26) in dilute asym¬ 
metric nuclear matter and on the structure of the cor¬ 
responding unstable modes. We start from an identity 
valid for any binary thermodynamical system 


/ / dpn \ _ (df^\ / dp.n \ _ 

\ dpp j j. \ dpn ) \ dpn } rp p^ \ dpp ) 

1 (( dl^p \ 

{l-y)p^ V dy j^p 


(29) 


where pq, pq {q = n,p) are respectively neutron/proton 
chemical potentials and densities. From our knowledge 
of the chemical potential on each isotherm, Eq.(9), we 
can easily compute the limits of the instability region in 
the T, pb plane for dilute asymmetric NM in the two 
choices, without and with the S meson. 

The results are reported in Fig.7. The 6 inclusion (solid 
lines) appears not much affecting the instability limits 
even at relatively large asymmetry a = 0.8 {N ~ 9Z). 
We notice just a small reduction and a shift to the left 
(lower densities) of the whole region: this can be un¬ 
derstood in terms of the larger symmetry repulsion, see 
Fig.3. 




for various asymmetries. Dashed lines: (p). Solid lines: (p+5) 
[for a 7 ^ 0]. 

We can understand the relatively small <5 effect on the 
stability border just remembering that for low densities, 
well below po; the symmetry term has roughly the same 
linear behaviour in both (p) and (p + <5) schemes, fixed by 
the 04 parameter (see the discussion after Eq.(20)). Only 
for very large asymmetries it appears relatively easier in 
the 6 case to be in the stable liquid phase. 

A rather larger difference can be seen in the behaviour 
of the quantity Eq.(26) inside the instability region. This 
is plotted in Fig.8 for various asymmetries at zero tem¬ 
perature and in Fig.9 at various temperatures for a fixed 
asymmetry a = 0.5 {N = 3Z). The solid curves (with 
(5-coupling) are systematically above the dashed ones, sig¬ 
nature of a weaker instability. 



FIG. 8. The quantity Eq.(26) inside the instability region 
at T = 0 and various asymmetries. Dashed lines: (p). Solid 
lines: {p + 5). 



FIG. 9. The same as in Fig.8 at various temperatures 
and fixed asymmetry a = 0.5. Dashed lines: (p). Solid lines: 
ip + S). 


In order to better understand the origin of such effect 
we study also the structure of the corresponding unsta¬ 
ble modes. We follow the Landau dispersion relation ap¬ 
proach to small amplitude oscillations in Fermi liquids 
31,^]. For a two component {n,p) matter the 
interaction is characterized by the Landau parameters 
Fq’'^ , {q,q') = (n,p) defined as 


N,iT) 


dpq 


dpq 


= Sa 


Fr 


(30) 


where Nq{T) represents the single particle level density at 
the Fermi energy. At zero temperature it has the simple 
form 


Nq 


kpqEp^ 


q = n,p. 


In the symmetric case {Fq^ = = Fq"), the 

Eqs.(27,28) correspond to the two Pomeranchuk instabil¬ 
ity conditions 


= F^'^ + F^P <-1 mechanical 
Fo“ = Fq"” - F^P < -I chemical. (31) 


^From the dispersion relations Fg will give the properties 
of density (isoscalar) modes and Fg of the concentration 
(isovector) modes. For asymmetric NM we have some 
corresponding generalized Landau parameters Fg^jF^^ 
which will characterize the new collective response (re¬ 
spectively isoscalarlike and isovectorlike). They can be 
expressed as fixed combinations of the Fg’^ for each 
baryon density, asymmetry and temperature. This trans¬ 
formation is reducing Eq.(29) to a ’’diagonal” form [p^ 


(l + FJg)(l + F“^) 


4 / y 

(1 - 2/)p2 \Nn + Np) 
\dp J rp y \ Oy J rp p 


( 32 ) 
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As already discussed in the unstable region of dilute 
asymmetric NAd we have isoscalarlike unstable modes 
and so 1 + < 0, while the combination 1 + Fq^ will 

stay positive. In Fig.10 we report the full calculation 
of these two quantities in the unstable region at zero 
temperature, for asymmetry a = 0.5, with and without 
the S coupling. 



inside the instability region at zero temperature and asym¬ 
metry a = 0.5 (N = 3Z). Dashed lines: (p). Solid lines: 
(P + S). 

The 6 meson is almost not affecting at all the unstable 
mode, given by the Fq^ parameter, and so the limits of 
the instability region, for a = 0.5, are not changed (see 
the Figs.7,8,9). We have a much larger effect on the Fq^ 
parameter which describe ’’stable” isovectorlike modes, 
that actually can propagate as good zero-sound collective 
motions since > 0. This can be expected from the 
isovector nature of the Fmeson. ^From here we get the 
main differences seen in Figs.8,9 for the quantity Eq.(26) 
inside the instability region, just using Eq.(32). 

Another interesting aspect of the comparison between 
Fig.8 and Fig.10 is the shift of the ’’maximum instability” 
density region. From the thermodynamical condition re¬ 
ported in Fig.8 it seems that the largest instability (the 
most negative value) is around ps = 0.06/m“^. In fact 
from Fig. 10 we see that the fastest unstable mode, cor¬ 
responding to the most negative Pomeranchuk condition 
for 1 -|- Fgg, is actually present for more dilute matter, 
around ps = 0.02/m“^. This shows the importance of 
the linear response analysis. 

Finally the fact that the 5-coupling is mostly affecting 
the isovectorlike modes is of great interest for possible 
effects on the isovector Giant Dipole Resonances stud¬ 
ied around normal density within the RMF approach in 
asymmetric systems. 


V. CONCLUSION AND OUTLOOK 


Aim of this work has been to select genuine new effects 
on properties of the asymmetric nuclear matter due to the 
coupling to the charged 5-scalar meson introduced in a 
minimal way in a phenomenological hadronic field theory. 
Simple analytical results allow to clarify the interplay 
between the p and 5 isovector meson contributions. 

With respect to the Equation of State, when the 5 
is included the symmetry energy keeps an almost linear 
baryon density repulsive behaviour at sub-nuclear densi¬ 
ties while it starts to move to a roughly parabolic trend 
above the saturation density po- Such effect comes di¬ 
rectly from the general relativistic property of the weak¬ 
ening of the attractive scalar interactions at high den¬ 
sity. It represents the equivalent in the isovector chan¬ 
nel (interplay of p and 5 contributions) of the saturation 
mechanism of the symmetric NM (interplay of uj and tr 
contributions). All that gives a further fundamental snp- 
port to the introduction of the 5-channel in the symmetry 
energy evaluation. 

Such ”5 mechanism” for the symmetry energy leads 
to a more repulsive EOS for pure neutron matter at 
baryon densities roughly above 2po- This is the region 
where transitions to different forms of nuclear matter are 
expected and so the result appears quite stimulating. 

We have shown in details that the new 5 contribu¬ 
tions are not negligible for the slope parameter around 
po (the symmetry pressure) and absolutely essential for 
the curvature parameter (symmetry incompressibility). 
The proton-neutron effective mass splitting is also di¬ 
rectly given by the 5 coupling and appears to be of the 
order of 15% at normal nuclear density (for a N = 3Z 
asymmetry). 

The possibility of an experimental observation of such 
effects is suggested. We list here some sensitive observ¬ 
ables: 


• Neutron distributions in n—rich nuclei (stable and 
unstable), ||25|-|^ 

• Assessment of drip-line stabilities 


• Bulk densities and incompressibility modulus in 
asymmetric nuclei [p^^2p^ ] 




Transport properties in Radioactive Beam colli¬ 
sions at intermediate energies: dissipative mech¬ 
anisms, fast nucleon emission, collective flows 


|^j2§3|||,||,|§ 


Effects on the new critical properties of warm ANAI, 
mixing of mechanical and chemical instabilities and 
isospin distillation, are also presented. The border of the 
instability regions as well as the nature of the unstable 
fluctuations in dilute asymmetric matter are not much 
affected. Indeed in the low density region the symmetry 
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energy has a very similar behaviour with and without 
the 6 meson. This reduces the possibility of the obser¬ 
vation of (5-coupling effects from measurements of isospin 
effects in multifragmentation events in asymmetric heavy 
ion collisions. We note however that other fragment pro¬ 
duction mechanisms, like ’’neck fragmentation”, are ex¬ 
pected to be very sensitive to the symmetry potential 
around normal density since clusters are now formed in 
regions in contact with the ’’spectator” matter 
We can predict some interesting effects on this kind of 
events, that actually represent a quite large part of the 
fragment production cross section. 

,^From the study of the 6 influence on the Landau 
parameters in the low baryon density region we note 
however an interesting effect on the collective response 
of asymmetric NM. While the unstable {isoscalarlike) 
modes are almost not modified, the isovectorlike ones, 
which appear to have a good stable zero-sound propaga¬ 
tion, are quite sensitive to the introduction of the 5. This 
suggests some importance of the 5 meson in the RMF 
description of isovector giant resonances in asymmetric 
nuclei. 

In conclusion we would like to add a comment on the 
limits of the Hartree scheme used in this work, as well as 
in most RMF calculations. Indeed it is well known that 
exchange terms can give contributions in the isovector 
scalar channel also in absence of explicit (5-couplings, see 
ref. and ref.s therein. The comparison with Dirac- 
Brueckner-Hartree-Fock results of ref. and some Non- 
Linear-Hartree-Fock calculations performed in |Q , , 

are showing that an effective S coupling parameter fs ~ 
2.5/m^, almost constant in a wide density range, is well 
accounting for the Fock term contributions. This makes 
us more confident on the reliability of the quantitative 
evaluation of the 6 effects discussed in this paper. 
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